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ABSTRACT 


^  A  new  method  of  estimating  wall  interference  in  unsteady  flow 
is  presented.  This  method  is  valid  for  subcritical  flow  within  the  accuracy  of 
the  linearized  small  disturbance  theory.  The  main  feature  of  the  method  is 
the  use  of  measured  pressure  along  lines  in  the  flow  direction  near  the  tunnel 
walls.  This  method  is  particularly  effective  in  a  tunnel  with  ventilated  walls 
because  it  does  not  require  the  representation  of  wall  characteristics  with 
unreliable  mathematical  expressions.  Results  of  some  numerical  examples 
indicate  that  the  new  method  produces  satisfactory  results  except  for  cases 
when  the  reduced  frequencies  are  dose  to  the  tunnel  resonance  frequencies. 

For  the  case  of  an  airfoil  in  pitching  motion,  the  method  has  been  used  to  , 
correct  the  amplitude  of  the  angle  of  attack  and  the  time  lag  in  the  motion.  A r~ 
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RESUME 

Une  nouveUe  methode  d ’estimation  de  la  perturbation  des 
ecoulements  instationnaires  par  des  parois  est  presentee.  Elle  s'applique  aux 
4coulements  subcritiques  dans  I’intervalle  de  precision  permis  par  la  thAorie 
des  petites  perturbations  lin6aris£es.  Sa  principale  caract&ristique  est  l’emploi 
de  mesures  de  la  pression  tangentielle  a  l’6coulement  pres  des  parois  de  la 
soufflerie.  La  m&hode  est  particuliArement  utile  dans  Ies  souffleries  a  parois 
ventilAes  car  elle  ne  n^cessite  aucune  representation  math£matique,  souvent 
erronee,  des  caractAristiques  des  parois.  Les  resultats  de  quelques  exemples 
numlriques  montrent  que  la  nouvelle  methode  donne  des  resultats  satis- 
faisants,  sauf  loraque  les  frequences  rMuites  appro chent  les  frequences  de 
r£sonnance  de  la  soufflerie.  Pour  un  profit  a£rodynamique  en  tangage,  la 
mAthode  a  permis  de  corriger  Tangle  d’attaque  et  le  retard  dans  le  motive- 
ment. 
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NOTATIONS 


Notation*  Definition 

a_  sound  speed  of  undisturbed  uniform  flow 

Cp  pressure  coefficient 

Cpo ,  Cpc ,  Cp,  defined  in  Equation  (22) 

Cpu  pressure  coefficient  at  z  =  H 


<&.<&•<£  Equation  (56) 
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Equation  (57) 
lift  coefficient 
defined  in  Equation  (92) 

|  CL  |  due  to  unit  pitching  motion 
defined  in  Equation  (27) 

defined  in  Equation  (21) 
defined  in  Equation  (53) 
airfoil  chord  length 
defined  in  Equation  (78) 

A 

-  0H 

tunnel  semi-height 

the  Struve  function 

amplitude  of  plunging  motion 

imaginary  number 

defined  in  Equation  (74)  or  (85) 

Bessel  function 

defined  in  Equation  (59) 

defined  in  Equation  (34) 

defined  in  Equations  (36)  and  (37) 
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Notations 
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Definition 
an  integer 

cd> 

reduced  frequency  *  — — 

zu^ 

a  large  positive  number 
uniform  flow  Mach  number 
an  integer 

normal  outward  co-ordinate  to  dS2  or  an  integer 

die  Newmann  function,  Equation  (A-21) 

defined  in  Equation  (A-10) 

a  large  positive  number 

a  large  positive  number 

defined  in  Equation  (80) 

defined  in  Equation  (83) 

mean  camber  of  airfoil 

transformed  time  variable,  Equation  (5) 

time 

uniform  flow  speed 
upwash  along  the  z-axis 
defined  in  Equation  (66) 
defined  in  Equation  (39) 
defined  in  Equation  (60) 
defined  in  Equations  (63)  and  (64) 
defined  in  Equation  (66) 
defined  in  Equation  (38) 
defined  in  Equation  (73) 

distance  measured  from  the  pitch  axis  in  the  uniform  flow  direction 
(vii) 
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Notations 
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*L 

A 

*T 

XL.*T 

A 
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a,  a 


“o.3» 


“a.  “a 
0 
7 

ACpc  i  ACps 

ACpc ,  ACpc ,  ACps 

ACp 

ACP 

AK 

AKj,  AK; 

AW 

A«a 

At 

d£2 

an, 

3S2r 
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Definition 

transformed  x,  Equation  (5) 

leading  edge  x-co-ordinate 

trailing  edge  x-co-ordinate 

transformed  xL  and  xT 

distance  upward  from  x-axis 

transformed  z,  Equation  (5) 

angle  of  attack 

time-averaged  angle  of  attack 

amplitude  of  angle  of  attack  of  pitching  motion 

defined  in  Equation  (4) 

Euler  constant  (=0.5772156649) 
defined  in  Equation  (53) 
defined  in  Equation  (32) 

pressure  coefficient  difference  between  the  upper  and  lower  airfoil  surfaces 

defined  in  Equation  (71) 

defined  in  Equation  (70) 

defined  in  Equations  (61)  and  (62) 

defined  in  Equation  (69) 

"  “A  "  “A 

time  lag  of  pitching  motion,  Equation  (1) 
boundary  of  SI 

a  small  circle  with  the  radius  p  and  the  center  at  (x,z) 
a  big  circle  with  the  radius  R  and  the  center  at  (x,z) 
a  small  positive  number 
integral  variables 
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Notations 

K 

(S.f) 
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Definition 

defined  in  Equation  (15) 

(x*) 

a  small  positive  number 
full  velocity  potential  functions 

small  perturbation  velocity  potential  functions,  Equations  (42)  and  (43) 
defined  in  Equation  (8) 
defined  in  Equation  (10) 
transformed  0,  Equation  (43) 
defined  in  Equation  (45) 
defined  in  Equation  (14) 
phase  lag  in  lift,  Equation  (91) 
analytical  real  valued  function  in  ft.  Equation  (47) 
defined  in  Equation  (17) 
integral  variable 
a  part  of  flow  field 

angular  velocity  of  pitching  or  plunging  motion 

transformed  do,  Equation  (9) 

the  Landau’s  symbol.  Equations  (A-30)  and  (B-9) 

„  1  for  x  >  0 
-1  forx<0 

averaged  value  over  the  airfoil 
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A  NEW  METHOD  OF  ESTIMATING  WIND  TUNNEL  WALL  INTERFERENCE 
IN  THE  UNSTEADY  TWO-DIMENSIONAL  FLOW 


1.  INTRODUCTION 

The  unsteady  wall  interference  problem  has  not  been  as  thoroughly  investigated  as  the 
steady  case  (Refs.  1  -  5).  A  new  method  of  estimating  the  wall  interference  in  unsteady  subsonic  flow 
is  presented  in  this  paper.  The  method  requires  the  time-dependent  pressure  along  lines  in  the  flow 
direction  near  the  tunnel  walls  and  on  the  airfoil  model  to  be  measured.  The  pressure  along  the  lines 
will  be  that  on  the  walls  if  the  tunnel  walls  are  solid.  This  method  is  particularly  effective  if  the  walls 
are  ventilated,  since  the  method  is  valid  no  matter  how  complicated  the  wall  characteristics  may  be. 
The  merit  of  this  approach  is  similar  to  that  for  steady  flow  which  were  developed  by  several  re¬ 
searchers  (Refs.  6  - 10).  There  is  only  a  slight  difference  between  the  steady  and  unsteady  wall  inter¬ 
ference  problem,  and  that  lies  in  the  fundamental  equation  for  the  small  perturbation  potential. 

In  incompressible  flow,  even  this  difference  vanishes  and  the  wall  corrections  can  be  obtained  using 
similar  method  developed  for  steady  flow.  The  new  method  will  first  be  presented  in  the  analytical 
way  and  then  some  examples  of  the  corrections  for  a  pitching  airfoil  will  be  calculated.  The  data 
used  here  are  not  obtained  by  tunnel  tests  but  by  numerical  means. 


2.  ANALYSIS 


Suppose  a  thin  airfoil  is  in  pitching  motion  with  a  constant  angular  velocity  in  free  air. 

In  this  case  it  is  convenient  to  locate  the  origin  of  a  space  co-ordinate  system  (x,z)  at  the  pitch  axis, 
where  x  is  the  distance  measured  from  this  axis  along  the  uniform  flow  direction  and  z  is  the  distance 
perpendicular  to  the  x  axis.  The  time  At  is  defined  as  the  time  when  the  time-dependent  averaged 
positions  of  the  airfoil  coincide.  The  angle  of  attack,  a,  (Fig.  1)  for  this  airfoil  can  be  written  as 
follows: 

<*(t)  =  <X0  +  0CA-sin.  [w(t  -a?)  |  .  (1 

In  order  to  simplify  the  analysis,  the  flow  is  assumed  to  be  inviscid,  irrotational  and 
subcritical  everywhere.  Furthermore  the  thin  airfoil  approximation  is  also  assumed  to  be  valid.  By  ^ 
the  aid  of  these  assumptions,  the  fundamental  equation  for  the  small  disturbance  potential  function  4> 
can  be  written  as 


(1-Mi)- 


+  -  2 


01 


A 


(2) 


where  <j>  is  defined  as  follows: 

4  ~  X  +  0  . 

v 

0(x,z)  is  the  full  velocity  potential  function.  and  a„  are  the  Mach  number  and  the  sonic  speed 
for  the  undisturbed  flow  infinite  upstream.  Variables  are  introduced  as  follows: 


=  c  u„-  ,  p  =  Vl-Mi  , 

ond  (x,  i  ,t  )  =  (  ox,  j.%  , 


where  c  is  the  chord  length  of  the  airfoil.  Equation  (2)  becomes 


(4) 


(5) 


-2- 


where  the  operator  V2  is 


v*  .  +j: 

dx 1  di7  ‘ 


(6) 


(7> 


Suppose  the  flow  disturbed  by  the  oscillating  airfoil  with  constant  angular  velocity  to  is 
also  periodical  in  time  with  the  same  angular  velocity  everywhere,  within  the  accuracy  of  the  small 

'S' 

disturbance  theory,  then  <t>  can  be  written 


(8) 


Defining  to  as 


gives 


CO 


(9) 


<p(.X,Z)  =  (J)g(cc,Z)  +  <ps(Z,Z)- sjn(cot)  +  <pc(.z,z)- COS(vt  )  J  (10) 


where  0O,  </>c  and  <p%  correspond  to  0Q,  <f>c  and  <j>%  in  Equation  (4). 
Substituting  Equation  (10)  for  0  into  Equation  (6)  gives 


v2$,  -  0  J 

(11) 

V2&  -  +  =0, 

(12) 

and 

V%  ♦  +"■•?.  =  o. 

(12a) 

Introducing  the  new  function  <b*  defined  as 

0  =  C0c  +  #£,  )•  e  .  (13) 


Equations  (12)  and  (12a)  can  be  combined  into  one  equation: 

->•  -  0  , 

*  •  <°/p  ■ 

From  Green ’s  theorem,  the  following  expression  for  can  be  obtained: 


(14) 


where 


(15) 


-3- 


(16) 


where  0  is  a  real  valued  analytical  function  throughout  the  flow  field  Cl  and  312  is  its  boundary  and 
n  is  the  outward  normal  co-ordinate  to  312.  The  following  form  for  i p  is  considered  here: 

VffejAtXjZ)  •  [  -=L=  •  sin{tcl^-z\{TT2  }-cos[KCt-x)e  )  de 

^  >o  v-e* 

VflM  6 


Cos{  <i^-x)G]de 


(17) 


is  continuous  and  differentiable  two  times  with  $  and  f  except  the  point  (x,z)  and  satisfies  the 
following  equation  except  at  that  point: 


vzij/  +  ty  =  0  . 


(18) 


Define  S 2  as  the  full  space  except  the  point  (x,z)  and  the  x  axis  downstream  from  the  leading  edge, 
Equation  (16)  reduces  to  the  following  expression: 

=  - ("[$«]!•  |f| 

1  -r  dE.  (£  =  /»  * 


where 


W-  =  (irn  {  0(4,4)-  0(4, -4)  j  . 


(19) 


(20) 


Appendix  A  can  be  referred  to  for  the  detailed  derivation  of  the  above  equation.  xL  is  the  point  which 
corresponds  to  the  leading  edge  x-co-ordinate  xL  in  Equation  (5). 


O'  A 

The  dynamic  pressure  coefficient  Cp(x,z,t)  can  be  calculated  from  0  as  follows: 


(21) 


Substituting  Equation  (8)  for  0  in  the  above  equation  gives 

A  A  .  ^  _  A  A  .  ^  K  A  A  A 


where 


and 


C,  A  A  N  _  .A  A  .  O'  .A  A  A  A  A  A  A  A 

p(x,  Z)  =  Cp0  c-x,  X  )  h  Cf,c(Zj*J'C0s(u>t)  +  CTs  (*,*)■  sin  (cot  )  , 


(22> 

(23) 

(24) 

(25) 
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With  the  aid  of  Equation  (14),  Equations  (24)  and  (25)  can  be  combined  into  one  expression: 

C-pCi,*, )  =  -  -4-'  C  ~  •  —  — •  i 

T  /i  ^  p>  3X  Mo,  '  V  • 


where 


Cp  —  CPr  +  i  Ci 


Solving  Equation  (26)  for  0*  with  the  boundary  condition: 


0*(x,z)  =  0  infinite  upstream. 


<p'(z,v--§  |  CpCcs„i,)-e 


From  Equation  (20)  follows 


where 


>*l!  =  -yf2[Vp(cO*e~i(^^~  x^d£, 

'  —  rv\ 


CP(i)  =  ACkCx )  +  i ACPS(*>  ■ 


ACpc  and  ACps  are  defined  as 


AC?(lb  =  ACj>0(x)  +  4  C?c(z) cc>s(u>t)  +  ACpJ>*)si*(u>t)  y  02) 

where  ACp(x,t)  is  the  time-dependent  pressure  coefficient  difference  at  x  between  the  upper  and  lower 
surfaces  of  the  airfoil. 

^  Differentiating  both  sides  of  Equation  (19)  with  respect  to  z  and  substituting  Equation  (30) 
for  [0*]  *  in  Equation  (19)  after  letting  z  tend  to  zero  give  the  following  expression 

=  y.p.j  CP(ce,)-R(e,-x)dt,  ,  *33) 

K(x)  «Rrw  +  iRiU)  .  ,34> 


and  k  is  the  reduced  frequency  defined  as 


*  =  £*  . 

2Uoo 


(35) 
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^  *># 

K,  and  K,  are  also  defined  as 

/V 


A 
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/V  jk  ^  A 

Kr(x)  =  f -cosUfc*)  -  *><**/ 2^  (^.-78*)*] 

““  A 


J^5  jj  jfj  -  0- ai’l2pr(e+MJ->.)<‘> 

4^y,  °  *”/2  fa-***!* 
-  % 


-  -&•  • »/. tt  J *)  ( 


7tyS 


— ^  ■ — ~—Jdn 

±-*r,uAt)-  £.1.00(2  2 ) 

|K.(  &  -0-<*njiy<e-fu*}<<t 


^  A 

KiCoc)  =  _  A. 

4Tl 

M,  n  m 

47TA 
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_  £ 


tu 

KJ1 


M, 

;  *■ 


)xf<t& 


w(x)  has  been  defined  as 


(36) 


(37) 


'*■'  A  A  •  A.  A 

W  (■*  )  =  Wt(x)  +  IW.ix), 


where  the  upwash  w(x,t)  along  the  airfoil  is  expressed  as  follows: 

W(x,t)  =  Wc  (X  )•  cos  (CO  t)  ■*■  Ws(a)-S»n(a?t )  . 

Because  the  airfoil  is  in  pitching  motion  expressed  in  Equation  (1),  wc  and  w,  can  be  calculated 
easily  and  given  as 


(38) 


(39) 


W(x)  =  -<XA  {2ft  ^  cos(co^)-sin(£de)jH0^cos(wd?)+2ft*  s»'n(&>d3)j.  (40) 


-6- 


From  Equation  (33)  the  pressure  distribution  on  the  same  airfoil  model  in  a  different  pitching  motion 
in  a  wind  tunnel  can  be  obtained.  The  uniform  flow  condition  infinite  upstream  in  the  tunnel  is  the 
same  as  in  the  free  air  and  the  angular  velocity  of  the  motion  and  the  fixed  point  are  the  same  and 
the  same  xz  co-ordinate  system  can  be  used.  However,  the  time  origin  is  not  defined  in  the  same  way 
as  in  the  free  air.  In  this  case,  the  time  when  the  position  of  the  airfoil  coincides  with  its  time-averaged 
position  is  defined  as  t  =  0.  The  angle  of  attack  a  for  the  airfoil  in  the  tunnel  can  be  written  as  follows: 

Q((t)  =  Of o  +-  QfA-s/n  (<£{.)  .  (41) 

The  half  height  of  the  tunnel  is  H  and  the  fixed  point  for  the  pitching  motion  is  on  the  tunnel  center 
line  in  this  paper  (see  Fig.  2).  The  small  disturbance  potential  function  0  for  the  flow  field  in  the 

tunnel  disturbed  by  the  airfoil  also  satisfies  Equation  (2)  like  0.  Here  0  is  defined  from  the  full 
velocity  potential  function  4>  as  follows: 


A  A  * 


$  =  Uoo*  x  +  0  , 

(42) 

Introduce  a  new  function  0  related  to  0  by  the  following: 

P>2i>  =  C-Uo-$  . 

(43) 

This  results  in  the  following: 

+  1C2  •  (p*  =  0  , 

(44) 

.  CoH+ 

where  <P*  =  (  <pc  +  i  <ps )  '  e  * 

(45) 

0C  and  0S  are  the  time-dependent  cosine  and  sine  terms  similar  to  0C  and  0S  in  Equation  (10).  Using 
Green’s  theorem,  the  following  expression  for  0*  can  be  obtained: 


where  0  is  an  analytical  real-valued  function  and  the  other  conditions  have  been  described  in  Equa¬ 
tion  (16).  Let  0  be  written  as 

/  //**f 

'  L  'HKW'-wWl  <47> 

then  this  function  is  continuous  and  differentiable  two  times  with  (  and  f  except  at  the  point  (x,z) 
within  the  following  band: 


l<  <  2H  . 


(48) 


Also,  kH  *  nit  ;  n*0,  tl.12,  ...  (49) 

since  the  principal  value  of  the  integral  must  exist.  0  is  satisfied  by  the  equation 


+  Kl  fljr  -  0 , 


(50) 


Hence,  Equation  (46)  becomes 


Appendix  B  can  be  referred  to  for  the  detail  derivation.  In  the  same  way  as  0*, 


Ml  =  -  -j  (  CP(cO- 

'-00 


^  0 


(52 


where  A 

CP(X)  =  ACt>c(x)  +  i25CPs(0. 

r  (53; 

ACpc  and  AC„S  are  the  time-dependent  cosine  and  sine  components  of  the  pressure  difference  between 
the  upper  and  lower  airfoil  surface  similar  to  A Cpc  and  ACps.  As  in  the  expression  for  [<p*]  t, 

-oo  /  y 

where 


CH(z)  . 

Cpc»  Cps.  ACpc  and  ACps  are  the  time-dependent  cosine  and  sine  components.  Noting  the  pressure 
difference  between  the  upper  and  lower  wind  tunnel  walls  with  CpH  (x,t)  and  using  a  Fourier  series 
expansion  t  gives 

CpH^t)  =  +  Cpc^->*  COS  (it.)  +•  Cps(X)*  si»  (u>0  (56) 


(54) 


(55) 


and  hence. 


CH  (x)  =  Cpc  (*>  *  i  CpS  (£ ) .  (57) 

Differentiating  the  both  side  of  Equation  (51)  with  z  and  substituting  Equations  (52)  and  (55)  for 
(0*J  +  and  (0*1  “H  into  Equation  (51)  and  letting  z  tend  to  zero  give 

XEp  =  v.p.  (  CP(c ,  (58) 

)-o» 

where 


(59) 


ind 


KU)  =  (Kr+iKi)  +  (4Kr+  \AKx  ) , 
WA(«)  *  +  *  W,(*X 


(60) 
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From  Equation  (47)  it  follows 

A 

AKr(r)  =  ■£{coth(2kH)-l}cos(2%x) 

A 

“  V'P-  (  •  Cot(2  {2  x(»iM „)j<# 

>  0  v+ft"  /  r 

"  V‘P'  ( '  -cot(2  ^Hltf)-*n{2&*x(f-AL)}<t9 

4nr  ’»  S  -  Jj-  Z3  Z3 


~^pvp  S{  ,  <•« 

A 

1K|  (*■)  “  “  -  I  }•  *"*  (2  fix) 

+  wrk'vp  f  'Z^-°*(2^Hfl:r)-c°H2&-x(B+M~)l<l9 
*KP  ***£.  ;  ^ 

'  Wp  Vf  \[  jff  cot(?  tf-M-M.)  1  <« 

r°°  r-r- 


'  4*fl  V?ll  jfx  C0t^  “'l*  tf-M-*-)  1 <« 

*  1  °  M#>  ’ 

fi  (62, 

C'  aa 

4’  e  2fth-  (  cCtK(2^fi;i-  |  }-<M(2$Z) 

~  jnpy'Filt  #7^  -“>^(2  )  vn[l  I  (t*M„)  j  dg 

-  Jy?  v P-  f  ^  •  ax^  fhP)-^{2^'x,l)-M.4^ 

-fe^l 

'  1  M"  (63) 


.  |M| 


I,  •=£ 


kWa 
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and 


WiU)  = 

+  J M  *^/c  CvTec(z^HYfp).a>T{z  *&?  *(*+&.)}& 

~4*pV?  f'p  jr^l/hp)-  ^  t 2  ^  x  (*~*~)  I dfi 

+  41t/5  ft  f^L^1  H  mi cJt  (2  ffq  )+!$■**  {2  |?"x  ($+M~)  J  do 

'  jnp'Pj  f^r  ^frfFi )*i)  u*{i (Huflt.  (64) 

Some  examples  of  AK  and  WA  are  shown  in  Figures  3  to  8.  w(x)  has  been  defined  as 


W(x)  =  tyc  (% )  +  l  U/s  (x)  , 
where  the  upwash  w(x,t)  along  the  airfoil  is  expressed  as  follows: 


(65) 


1 


3F(x,t)==  UJo  (x)  +  WiCZ)'Cos(a>t)  +  Vs (."*)  •  s'n  (At)  .  (66) 

From  Equation  (41 )  it  follows  that 

W(x)  =  -Z%\  -iofo  _  (67) 


Subtracting  Equation  (33)  from  Equation  (67)  gives 

i  o.  (Zt  r 

^Sto(x)  -  ]  CP(c£,)-4K(£.-^<<£,-  CH(p£,)-H'A(£.-i)dE. 

ZL  '-oo 


+  v.p.  r,4CPCcE,)-K<E.-%M£,  , 

(68) 

where 

A  w(x)  =  w( x )  -  >v (^c )  , 

(69) 

=  AKr(*)  +  iAK  i_, 

(70) 

and 

ACP(x)  =  CP(*>  -  CP(x). 

(71) 

It  should  be  noted  here  that  CP  (x)  and  CH  (x)  and  w  (x)  can  be  measured  in  a  wind  tunnel  test  but 
£p(x)  and  aA  and  Ar  are  unknown.  The  unknown  quantities  in  Equation  (68)  are  ACP(x),  aA  with 
A r  being  an  arbitrary  parameter.  The  aim  in  the  wind  tunnel  wall  correction  is  to  obtain  interference- 
free  data  from  the  test  data.  Therefore,  it  is  better  for  CP  and  CP  to  be  dose  to  each  other  because 


there  is  the  possibility  that  the  difference  between  CP  and  CP  can  generate  completely  different  flow 
field  from  each  other  due  to  the  coupling  between  viscous  and  in  viscid  flow  regions  such  as  the 
separated  flow.  As  a  result,  aA  and  At  must  be  chosen  so  that  ACP  may  be  as  small  as  possible.  ACP  is 
a  function  of  z  and  the  magnitude  of  it  depends  on  its  definition.  Equation  (68)  reduces  to 

V-P*  f  ACP(c£,)‘Ra,-z)d£  =  +  ,  (72) 

'  X.  Us»  U  •• 


where 


»i(l) 

Uoo 


CP(c£)-4K«-*)^  +  [  CH(c^)  WAC^-x)dBi. 

r,  -oo 


In  this  paper,  aA  and  A r  are  chosen  so  that  the  following  function  may  be  a  minimum: 


i(5A,4*) 


/  *r 

1  j  IVf  +4H'|2<lt>  • 


(73) 


(74) 


It  will  be  studied  in  the  next  section  whether  these  values  of  aA  and  Ar  are  satisfactory.  Then 


ai 


(75) 


and 


ai_ 


(76) 


Prom  Equations  (40),  (67)  and  (69)  it  follows 


where 


tan  CcoA^) 


ffc  _ 

0(A(  t+4^-xl)  “  (  f s  +  2  ) 


(77) 


+ c  +  if 


*  > 


(78) 


and  (  )  means  the  averaged  value  over  the  airfoil  chord. 


Also, 

S*  -  tfA  =  + 


- - - k-ft-2fe-*fs 

1  +  4^-X* 


-  C  f,  t  )-cc>s(&A'?)]  . 


(79) 
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Because  the  values  of  Ar  and  aA  can  be  calculated  with  the  aid  of  Equations  (77)  and  (79),  the 
right-hand  side  of  Equation  (72)  can  be  estimated.  If  the  value  of  the  term  in  Equation  (72)  is  small 
enough  over  the  airfoil  to  be  neglected  within  the  accuracy  of  the  test  condition,  ACP  can  be  also 
neglected.  In  this  case,  the  tunnel  test  data  can  be  corrected  in  the  sense  of  wind  tunnel  wall  inter¬ 
ference.  However,  not  all  tunnel  test  data  can  be  corrected.  In  such  circumstances,  it  is  necessary  to 
solve  the  integral  equation  Equation  (72)  for  ACP(x)  and  to  get  CP(x)  with  the  aid  of  Equation  (71). 
However  such  additional  correction  is  not  recommended  because  it  is  desirable  that  tunnel  test  data 
are  affected  as  little  as  possible.  In  addition  to  the  correction  mensioned  in  this  paper,  the  correction 
corresponding  to  the  time-averaged  flow  must  be  made.  This  correction  is  the  same  as  that  for  steady 
flow.  The  terms  fc  and  fj  on  the  right-hand  side  of  Equation  (78)  can  be  calculated  from  Equation  (73) 
with  the  aid  of  Equations  (61)  to  (64). 

The  wall  interference  on  an  airfoil  in  plunging  motion  can  be  estimated  in  a  similar  way  as 
that  for  pitching  motion.  The  time  and  space  co-ordinate  system  is  the  same  as  in  th /^previous  case. 
The  z -co-ordinate  of  the  mean  camber  of  the  airfoil  model  in  the  motion  in  free  air,  Sc(x,t),  is 

Sc(x,t)  =  SCo(X)  +  V  sin(cv-(t  -d?  >)  .>  (80) 

where  SC()  means  the  mean  camber  of  the  airfoil  itself  and  hA  is  the  amplitude  of  the  plunging 
motion.  In  this  case  the  upward  velocity  at  a  point  on  the  airfoil  is 

V(x,t)  =  fiAU)*  cos{  co 
Using  this  definition  for  w(x),  it  follows 

VV  ("x )  —  •  oj  •  cos  CtoA’i )  +  I  •  u>  •  sin  (tod?  )  .  (82) 

Consider  the  same  airfoil  model  in  the  plunging  motion  installed  in  a  wind  tunnel.  The  angular  velocity 
is  also  the  same  as  to,  and  the  z-co-ordinate  of  the  mean  camber,  Sc(x,t),  is 

Sc(£,t  )  =  SC„(X)  +  flA-sin(wt).  (83) 

From  Equations  (65)  and  (66), 

VKx)  =  .  (84) 


By  the  same  method  used  in  deriving  corrections  for  pitching  motion,  the  corrections  for  the  plunging 
motion  can  be  estimated.  In  this  case,  the  estimation  function  I  in  Equation  (74)  becomes  a  function 
of  hA  and  Ar,  that  is, 

✓V  fXj  2 

=  I  W\tc£.)-t-4*r(cg>)l  dfc,  .  (85) 

This  function  must  be  minimum  with  respect  to  both  hA  and  Ar: 


and 


a  i.  n 

i*  "  0  ’  ,86) 

II  =  0. 

3A?  (87) 


A 
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From  Equations  (85)  to  (87), 


)  —  — - -  ,  (88) 

K+Zl-ll a/c 

and 

ft*  —  =  l  C0S(^4^)-J  j-8*  +  ic  C0S(u>A'l)  +  fs-sinO&A'fc)  ]  •  (89> 

3.  SOME  EXAMPLES  OF  CORRECTION 

Some  examples  of  wall  corrections  using  the  new  method  are  described  here.  It  should  be 
noticed  here  that  the  examples  indicate  only  the  process  of  using  the  method  and  its  limitations, 
but  that  they  cannot  prove  the  validity  of  the  method.  This  method  has  been  already  proved  analyti¬ 
cally  on  the  base  of  some  of  the  assumptions  mentioned  in  the  previous  section.  Many  experiments 
have  indicated  these  assumptions  are  reasonable  and  they  have  been  adopted  in  many  other  papers. 

For  simplicity  in  the  calculations,  the  wall  interference  in  an  open  tunnel  is  investigated. 

In  order  to  confirm  the  procedure  of  solving  Equation  (33)  to  be  correct,  the  pressure  distributions 
on  the  oscillating  airfoil  in  pitching  motion  are  calculated  using  a  numerical  method.  The  airfoil  is 
installed  in  an  open  tunnel  having  a  half-height  of  5.0.  The  uniform  flow  Mach  number  is  0.866. 
The  airfoil  chord  length  is  1.0  and  the  pitch  axis  is  at  the  mid-chord  point.  The  amplitude  of  the 
oscillating  angle  of  attack  is  1.0  degree.  Ten  values  of  the  reduced  frequency  are  adopted  ranging  from 
0.02  to  0.18  in  steps  of  0.02  plus  0.17  because  0.182  is  approximately  the  tunnel  resonance  frequency 
which  can  be  calculated  from  Equation  (49).  The  tunnel  resonance  frequency  for  the  open  tunnel 
cannot  satisfy  the  condition  given  in  Equation  (49).  The  amplitude  of  CL  due  to  unit  pitching  motion 
with  reduced  frequency  is  shown  in  Figure  9,  and  the  phase  lag  in  CL  is  shown  in  Figure  10.  These 
quantities  are  calculated  from  the  expression 

JClU  =  +  ,  (90) 

and 

<J>ld  =  Tanl(  Cls/Clc  )  >  (91) 

where 

CL(t)  -■  CLo  +-  Clc  C°sCwI)  1-  CLs- sin^t)  .  (92) 

The  two  figures,  Figures  9  and  10,  show  good  agreement  with  the  figures  in  Reference  3.  This  fact 
means  the  present  solving  process  is  right.  Then  suppose  the  data  obtained  by  this  numerical  way  to 
be  data  which  should  be  corrected  to  data  in  the  free  air.  In  this  case,  the  pressure  coefficient  on  the 
tunnel  walls  is  always  zero.  With  the  aid  of  Equations  (77)  and  (79),  the  incidence  amplitude  and  the 
phase  lag  in  free  air  corresponding  to  the  tunnel  flow  can  be  calculated  from  the  pressure  coefficient 
distribution  on  the  airfoil  in  the  open  tunnel.  Results  are  shown  in  Figures  11  and  12.  The  pressure 
coefficient  difference  between  the  upper  and  lower  surfaces  of  the  airfoil  is 

ACp(x,t)=  ddf>9(z)  +  t  4Cp*(J)- l)  ,  (93) 

which  has  been  defined  at  Equation  (53).  The  time-dependent  angle  of  attack  is  expressed: 

0(<t)  =  0.0/74533  •  s'»(£i  )  . 


(94) 
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From  Equation  (1)  follows 


a*  A 

oUt) 


<yA  ■  sin(cv-Ct  -4"5)/ 


(95) 


For  the  pressure  coefficient  difference  on  the  airfoil  in  the  free  air.  This  equation  can  be  approxi¬ 
mated  by  Equation  (93)  when  the  motion  of  the  oscillating  airfoil  is  defined  by  Equation  (95),  where 
aA  and  coAr  can  be  calculated  using  Equations  (77)  and  (79).  If  the  airfoil  in  the  free  air  oscillates 
in  the  same  way  as  in  the  tunnel,  that  is, 

ou!;  =  5 A  •  s/n(u>t)  J  (96) 

A 

the  time  variable  t  is  transformed  by 

A  7T 

t  =  £  +  AQ  .  (97) 

The  resultant  expression  for  ACp(x,t)  can  be  written  as: 

^  A  ^  /V. 

ACj,(z,  t  )  =  j  hCycjA  )■  cosCiaAZ)  t  ACps  (*)■  si"  (ujA?)}  •  cos(&  t  ) 

+  { -  ACfc(*)- v'n(u)  t  )  (98) 

In  this  way,  the  in-phase  and  out-of-phase  components  of  the  pressure  coefficient  difference  as^a 
function  of  the  oscillating  angle  of  attackcan  be  calculated.  By  solving  Equation  (33)  directly,  ACp 
can  be  obtained.  Comparison  between  ACp  obtained  from  Equation  (98)  and  Equation  (33)  provides 
a  simple  way  for  correcting  aA  and  At.  The  results  are  shown  in  Figures  13  to  14.  These  figures  show 
that  the  method  is  suitable  for  low  reduced  frequencies  but  fails  near  the  tunnel  resonance  frequen¬ 
cies.  The  time-dependent  lift  coefficient  and  pitching-moment  coefficient  about  the  pitch  axis  with 
the  angle  of  attack  are  also  shown  in  Figures  15  to  18.  This  example  gives  the  extreme  case  because 
the  uniform  flow  Mach  number  is  0.866.  The  next  example  is  the  same  as  in  the  previous  case  except 
for  the  Mach  number  and  the  tunnel  height.  The  uniform  flow  Mach  number  is  0.600  and  the  tunnel 
semi-height  is  4.0.  In  this  case,  the  tunnel  resonance  reduced  frequency  is  0.52360.  The  reduced 
frequency  was  varied  from  0.05  to  0.50  in  steps  of  0.05  and  0.52.  Plots  of  |  CL  |  and  <pLa  with  k 
are  shown  in  Figures  19  and  20.  The  corrections,  (aA  -  aA  )  and  wAr,  are  shown  in  Figures  21  and  22 
while  Figures  23  to  24  show  ACp  and  Figures  25  to  28  show  the  CL  vs  a  and  Cm  vs  a.  These  figures 
show  also  the  same  behaviour  as  in  the  previous  case. 


4.  CONCLUSIONS 

A  new  method  of  estimating  the  wind  tunnel  wall  interference  on  an  oscillating  airfoil  is 
presented.  Instead  of  expressing  the  wall  condition  in  unreliable  ways,  the  time-dependent  pressure 
distributions  on  the  oscillating  airfoil  and  near  the  tunnel  walls  in  the  flow  direction  are  measured. 

With  the  aid  of  the  measured  pressure  distributions,  the  corrections  to  the  incidence  amplitude  and  the 
phase  lag,  (uA r),  can  be  calculated  by  this  new  method.  This  method  is  very  effective  when  the  test  is 
performed  in  a  tunnel  with  ventilated  walls  because  the  unreliable  expressions  for  the  wall  character¬ 
istics  are  not  used  at  all  in  this  method.  The  corrections  are  satisfactory  except  at  reduced  frequency 
near  the  tunnel  resonance  frequency.  To  determine  whether  corrections  can  be  carried  on  by  this 
method,  it  is  necessary  to  check  the  induced  upwash  distribution  on  an  airfoil  due  to  the  walls. 
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FIG.  2:  AN  AIRFOIL  IN  PITCHING  MOTION  IN  A  TUNNEL 


-17  - 


WEIGHTING  FUNCTION  AK 


UNIFORM  MACH  NUMBER  =  0.60 
REDUCED  FREQUENCY  =  0.100 
TUNNEL  SEMI-HEIGHT  =  4.00c 


AKr 

AKi 


FUNCTION  WA 


G.B 


-23- 


Ll  FT  COEFFICIENT  DUE  TO  UNIT  PITCHING 

MOTION  WITH  REDUCED  FREQUENCY 
UNIFORM  FLOW  MACH  NUMBER  =  0.866 
CENTER  OF  PITCHING  MOTION  :  0.50 


A 


FIG.  9:  I  CL  |  vs  k 

<k 


-24- 


PHASE  LAG  IN  LIFT  DUE  TO  UNIT  PITCHING 

MOTION  WITH  REDUCED  FREQUENCY 
UNIFORM  FLOW  MACH  NUMBER  =  0.866 
o  CENTER  OF  PITCHING  MOTION  :  0.50 

o 

s- 

!  —  free  % 

o  OPEN  TUNNEL 

It 

SEMI-HEIGHT  =  5.00c 


O'  t  o 

oj 

fO 

I  • 


o  1 
J 


FIG.  10:  |0L|o  v*  k 


INCIDENCE  AMPLITUDE  INCREMENT  DUE  TO 

WALL  INTERFERNCE  WITH  REDUCED  FREQUENCY 
UNIFORM  FLOW  MACH  NUMBER  =  0.866 
CENTER  OF  PITCHING  MOTION  :  0.50 

o 


OPEN  TUNNEL 
SEMI-HEIGHT  =  5.00c 
INCIDENCE  AMPLITUDE  =  1.0° 


-26- 
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APPENDIX  A 

DERIVATION  OF  EQUATION  (19) 


Firstly  it  will  be  proved  that  *,z)  in  Equation  (17)  satisfies  the  Helmoltz  equation 
Equation  (18)  everywhere  except  at  a  point  (x,z).  Equation  (19)  will  be  obtained  later. 

i^(£,f;x,z)  in  Equation  (17)  is  bounded  everywhere  except  at  a  point  (x,z).  Unless  f  =  z, 
the  derivative  of  with  £  or  f  can  be  obtained  by  the  way  of  differentiating  the  two  integrands  on 
the  right-hand  side  of  Equation  (17)  directly  with  £  or  f .  So  i p  satisfies  Equation  (18)  unless  f  =  z. 
From  Equation  (17)  it  follows: 

=  ~sgr\(4-2?)-{  (  Cosl<\£,-Z\)l\Ip  }-Cos|K(£,-r)0  J  d0 

(A-l) 


•  Cjos[Ki(t>-z)&}cl8  I  ■ 


Equation  (A-l)  can  be  arranged  to  give: 


Vk\ 


K 

m 


■sr 


cos{  <(k~  x)6}  dd  tie*1*’  cos  x)&}  dp 


'  (? 

_  ]IP^i  j'  Msf  /c(4 -x)  9] ]  • 


(A-2) 


Because  the  first  two  integrals  at  the  right-hand  side  of  Equation  (A-2)  can  be  estimated  analytically 
easily,  Equation  (A-2)  reduces  to 


3y/ 

iZ 


I 


_  J_  \<-zl 


tt,- */*«-*/ 


•  COS  {  Kt$  -  X  )  j 


I  I  -cos(icJ4-^|y/_^*; )-cos{K(£>-z)fd0 


So 


Him  =  Him  =(?;£,  *  X  . 


( A-2a) 


(A-3) 


Because  \p  is  continuous  with  f  except  at  a  point  (x,z)  and  because  of  the  relation  (A-3),  the  derivative 
of  i/'  with  f  on  the  line  f  =  z  exists  and  its  value  is  also  0  unless  £  1  x.  The  derivative  of  with  f  is 
also  continuous  with  respect  to  f  except  at  a  point  (x,z).  The  second  derivative  of  \p  with  f  is  obtained 
from  Equation  (18): 


v  v 


£  ( I  _ 

27T  |  U-0*  •  s"*  I  *l< -tiihP  ]•  COii  ic(i-xj6  Jd& 
0 


at  e 


-KIC.-ZI 


cos  /  «(£  -x)9  }  <19  . 


(A-4) 
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From  Equation  (A-4),  the  second  derivative  is  symmetric  with  respect  to  (f-z).  Resultantly  the 
second  derivative  does  exist  and  its  value  is 


t  X 


a,..  =  tim 

because  of  the  same  reason  as  the  first  derivative  is  estimated.  From  Equation  (18),  follows 


(A-5) 


$ {  N'UcIZ-xl  )  ,  (A-6) 

where  N0  is  the  Newmann  function.  So  the  first  and  second  derivatives  with  £  exist  unless  £  =  x. 
Because  the  second  derivative  of  i p  with  £  at  f  *  z  is  symmetric  with  respect  to  (£-  z), 

vw  _  #  rl\ 

U*’  (A-7) 

It  is  of  course  that  the  second  derivatives  of  \p  with  both  £  and  f  are  bounded  even  if  f  tends  to  z 
unless  £  =  x.  Resultantly  i^(£,fpc,z)  can  satisfy  Equation  (18)  everywhere  except  for  a  point  (x,z). 


Now  suppose  a  small  circle  in  the  considered  space  (£,f )  which  has  a  radius  p  and  its  center 
at  (x,z).  An  arbitrary  point  (£,f)  on  the  circle  can  be  expressed: 


£,  =  X  -t  P  ■  CDs  it 

'  ) 

<  =  ?+/>•  sjnih  ■  f>  >  0 

If  n  is  the  inward  normal  co-ordinate  to  the  circle, 

d  ys  _  djjr 

dn  ~  df> 


(AS) 


(A-9) 


Therefore 


-  -  ft I '  t'cfl™*l0r'.osil'eo*s*>eidi> 


\  W-i 

where  P  is  bounded  for  both  p  and  This  circle  is  notated  with  3f2p. 

■2n. 


4  |f  P*ds  =  -p*(z„z,)-r [  ^ dv> 
hnP  dn  >0  J 


where  (x0 ,  z0 )  is  in  3£2p. 


Because 


p-0  '  )0  9f 


,2V 


tim  p  P(  P,^»  )  di(,  —  0  t 


(A-10) 


(All) 


(A-12) 


(A-13) 
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.  1 2lt  /  o® 

fc)  P  [^Is,n^1  I  (A-14) 

'  0  I  I  1  In 


and 


Him 

P-o 


( 27f  ,  c*  _  /  27T 

P  lp  S(  fg^' e  Kp|SlnKl1^' s,M|lcP(cL0S«:>*/d9]d^'=  i^rj  C03’v>ci'*  ■  (A'15) 


Then 


cx>s{ Kf>(a>sik)Q]ci&  +  CL(l°jit)  ,  (A-17) 


Jim  i  ^ •  <p*  ds  =  -<p*Cz,Z)  - 

P-o  >*1,  Jn 

Similarly  ip  on  the  circle  is 

_  J_  (°°  I  ^i^im 

W  2fC  Jf  IFi  e 

where  Q(p,i// )  is  bounded  for  both  p  and  \p.  Then 

1 L;  -  £<s  1  s  pM  1 0-i- 

where  M  is  the  maximum  value  of  d<p*/dn  on  d£lp.  From  Equation  (A-17)  follows 

'*1  i  *  lei. 

The  integral  at  the  right-hand  side  of  Equation  (A-19)  reduces  to 

JL I  (  _J_  I 

2Tf  J,  tfTTT  e  oos(KP(OQjv)BidB  I 


<  ±r » 

2T^  Id 


-pKlSl'n>|  3 


(A-16) 


(A-18) 


(A-19) 


/o^  ^  -d£  =  ^{He(»tplswltl)-/V6^/>/»>’«f/;]  2  (A-20) 

where  H0  is  the  Struve  function  and  bounded  for  any  value  of  its  variable.  The  Newmann  function  N0 


has  the  following  property : 


N,w  =  4 -J.<w 


2  V  H)  ,%?  ?  ± 

11 L  an'  (z)  k,m  J 


(A-21) 


where  J0  is  the  Bassel  function  and  bounded  for  x  and  y  is  the  Euler  constant  number,  (y  -  0.57721...) 
Then 


Km  p  ■  |  {j;  |  =  0  . 
P-0 


(A-22) 


Resultantly  from  Equation  (A-18)  follows 

1  s£*dsl-~°  ;  t  —  °  .  'a-2 

With  the  aid  of  the  thin  airfoil  approximation  Equation  (16)  becomes 

■?*«.<>= -/"[♦•■£-*.*£•]:*£  »•)«!. , 

where  3£2R  is  a  circle  with  the  radius  R  and  its  center  at  (x,z).  As  previously  mentioned,  ^  and  its 
first  derivative  with  f  are  continuous  except  for  a  point  (x,z).  So  the  first  integral  becomes 

30*/3f  means  the  upward  velocity  component  and  should  be  continuous  throughout  the  flow  field. 
But  it  must  coincide  with  the  vertical  component  of  the  airfoil  surface  velocity  from  xL  to  xT .  xL 
and  xT  are  the  transformed  leading  and  trailing  edge  co-ordinates,  respectively.  Because  the  airfoil 
model  is  oscillating  in  one  piece,  there  is  no  difference  in  the  vertical  velocity  between  the  upper  and 
lower  surfaces.  Then 

[  ]-  =  0  ;  X  £  [xL>oo)  .  (A-2 

As  a  result  the  first  integral  of  Equation  (A-24)  becomes 


From  Equations  (A-10)  and  (A-17)  it  follows  that 


I  )  ’  C R  -*°°  3  . 


1  Vlf-sl  -  0(3=0  ;  C  R  -~]  , 


where  ‘O’  is  the  Landau’s  symbol  and  (A-28)  means 


iim  Vr  .1  m  I  =  c 
R -oo  m  1  9PW  C  • 


l  ^Ip-r  -  o(  j 


*?\  I  =0(1.) 

*P  Ip-Rl  UKfRJ' 


R-oo  Jao„  drt  dn  '  ' 
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In  this  case  Equation  (A-24)  becomes  the  same  equation  as  Equation  (19).  However  CP(x)  obtained 
from  the  integral  equation  Equation  (33)  based  on  Equation  (19)  is  limited  because  0*  corresponding 
to  dP  always  satisfies  the  conditions  (A-31  and  (A-32).  Experimentally  such  CP  as  is  obtained  from 
Equation  (33)  seems  a  good  approximation  for  the  flow  considered.  Both  conditions  (A-31)  and 
(A-32)  are  assumed  to  be  satisfied  in  this  paper. 
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APPENDIX  B 

DERIVATION  OF  EQUATION  (51) 


The  integral  function  on  the  right-hand  side  of  Equation  (47)  is  analytical  everywhere 
including  the  point  (x,z).  The  function  also  satisfy  the  Hermholtz  equation  Equation  (50).  As  a  result 
4/  is  also  a  solution  of  Equation  (50).  Because  \p  can  be  expressed  in  the  form  of  the  sum  of  \p  and  an 
analytical  function  everywhere,  Equation  (46)  reduces  to  the  following  form  similarly  to  Equa¬ 
tion  (16):  (see  Fig.  29) 

=  -jVf 

From  the  same  reasoning  as  mentioned  in  Appendix  A, 

C  H  Jt.  =  0  »  x  e  [Xl.oo)  . 


(B-l) 


The  last  term  in  Equation  (B-l)  vanishes.  Infinitely  upstream  there  is  no  disturbance  so  that 


h.  L(*\f -*-$f  -o 

From  Equation  (47)  it  follows 

=  ;  cr 


where  ‘O’  is  the  Landau’s  symbol  in  Equation  (B-4)  and 


oo 
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tlm  VR-I  01  I  =  C 

R~<*>  '  1  t-R  I  > 


(B-2) 


(B-3) 


(B-4) 


(B-5) 


where  c  is  a  constant  number.  Similarly 


Let 


1  | 

X 

0 

;  C  R  — «» J . 

(B-6) 

-  o(fR) 

;  hr— , 

(B-7) 

I  1 

1  a*1*.*! 

-  o({r) 

;  Lr  —  -0  , 

(B-8) 

and 
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where  ‘O’  is  also  the  Landau’s  symbol  and  Equation  (B-7)  becomes 


£im.  1 1  _ 


(B-9) 


Then 


H 


fR 


*  0 


-H 


oo 


(B-10) 


Experimentally,  the  conditions  (B-7)  and  (B-8)  are  always  satisfied  except  for  the  case  of  tunnel 
resonance.  Equation  (B-l)  then  becomes  Equation  (51). 
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